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We present constraints on the number of relativistic species from a joint analysis of cosmic mi- 
crowave background (CMB) fluctuations and light element abundances (helium and deuterium) 
compared to big bang nucleosynthesis (BBN) predictions. Our BBN calculations include updates of 
nuclear rates in light of recent experimental and theoretical information, with the most significant 
change occuring for the d(p, 7)^He cross section. We calculate a likelihood function for BBN theory 
and observations that accounts for both observational errors and nuclear rate uncertainties and can 
be easily embedded in cosmological parameter fitting. We then demonstrate that CMB and BBN 
are in good agreement, suggesting that the number of relativistic species did not change between 
the time of BBN and the time of recombination. The level of agreement between BBN and CMB, as 
well as the agreement with the standard model of particle physics, depends somewhat on systematic 
differences among determinations of the primordial helium abundance. We demonstrate that in- 
teresting constraints can be derived combining only CMB and D/H observations with BBN theory, 
suggesting that an improved D/H constraint would be an extremely valuable probe of cosmology. 

PACS numbers: 98.80.-k, 26.35.-l-c, 98.70.Vc, 14.60.St, 25.10.+S 



I. INTRODUCTION 

The sensitivity of early-universe cosmology to the num- 
ber and properties of light particle species was recognized 
early in the development of the hot big bang model 

The expansion rate of the universe at early times 
increases with number of relativistic particle species in 
thermal equilibrium, and this in turn sets timescales for 
big bang nucleosynthesis (BBN). One can then use the 
BBN yields of light nuclei to constrain the number of light 
species quantitatively Q. Constraints on the expansion 
rate are therefore often expressed in units of an effective 
number of light fermion species A^cff, so that for the stan- 
dard model of particle physics A^off ~ 3, while iVcft ^ 3 
corresponds to some altered expansion rate that might be 
caused by additional neutrino species or by some other 
effect (like a nonzero neutrino chemical potential or vi- 
olations of general relativity) [t^ . Recent BBN-derived 
constraints may be found in Refs. [^-Q. The cosmic mi- 
crowave background (CMB) is sensitive to the expansion 
rate much later: at the times of equal matter and radi- 
ation densities and of hydrogen recombination. Recent 
measurements of the CMB power spectrum now allow 
interesting constraints on A'eff at those times (ol-fl^. 

Constraints on light species from BBN and the CMB 
have been examined many times as both observational 
data and theoretical calculations have improved. The 
consistency or not of these two types of constraints is an 
important question, particularly since they probe expan- 
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sion rates at widely separated times in the evolution of 
the universe. In the present work, we check for consis- 
tency between the two data sets, examine the sensitivity 
of that consistency to model assumptions, and then at- 
tempt to place the strongest constraint possible on iVefi 
by combining the two types of data. 

The BBN yields - especially of deuterium - depend not 
just on iVcff but also on the cosmological baryon density, 
which in units of the critical density of the universe is 
rif, = SnGpB/Sh^. {h is the Hubble constant in units 
of 100 kms^^ Mpc~^, G is Newton's constant, and pg is 
the baryon density expressed as mass per unit volume.) 
Prior to the precise inference of fib/i^ from CMB mea- 
surements, the strongest constraint on N^g came from 
using the deuterium-to- hydrogen number ratio D/H to 
restrict flth^ and then exploiting the TVcff dependence of 
the primordial helium mass fraction Yp. However, D/H 
has its own dependence on A^^eff [1, 113] ; a strong external 
constraint on flbh^ allows BBN constraints on A^off that 
are independent of Yp. We show below that, particularly 
because precise measurements of Yp are difficult, the con- 
straint on iVoff from D/H is competitive with that from 

In this paper, we focus on three issues: a full update 
of the BBN reaction rates including a revision to the 
(i(p, 7)'^He rate based on recent theoretical calculations; 
an implementation of BBN constraints that includes both 
observational and theoretical uncertainties within Cos- 
moMC [l3| (a popular Markov chain Monte Carlo fitting 
program used for cosmological data); and the utility of 
D/H as a probe of new physics that is independent of 
uncertainties in the primordial helium abundance. 
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II. BIG-BANG NUCLEOSYNTHESIS AND N^, 



The reaction p{n, 'y)d 



BBN constraints on physics and cosmology require ac- 
curacy in both computed BBN yields and observed abun- 
dances. With only standard-model physics {N^=3, no 
late-decaying particles, negligible neutrino chemical po- 
tentials, etc.), BBN is a one-parameter model, depending 
on the baryon density. Since Qbh"^ = 0.02249(55) [9| is 
now known to ^ 3% from the CMB alone, this "stan- 
dard" BBN is a high-precision theory with no free pa- 
rameters. The predicted yp has a precision of ~ 0.2%, 
while yields of D, ^He, and ^Li are predicted with respec- 
tive precisions of roughly 5%, 4%, and 8% (computed 
below and e.g. in Refs. |15l - [l9| ). To constrain devia- 
tions from this standard model, we need both accurate 
BBN calculations based on up-to-date nuclear inputs and 
observationally-inf erred primordial abundances. 



A. The BBN calculation 

We compute BBN light-nuclide yields using a modi- 
fied version of the Wagoner/Kawano code described in 
Ref. 20]. We adjusted the timestep controls to com- 
pute yields more precisely than 0.1% for given nuclear 
rates, and we used the Monte Carlo techniques described 
in Ref. [13] to propagate nuclear cross section data into 
yields and yield uncertainties. For the centroid of our Yp 
predictions we used the independent BBN code of Lopez 
and Turner [2l'|. It accounts for several small effects in 
weak-interaction physics and the equation of state and 
achieves a theoretical precision (for a given neutron mean 
lifetime t„) of AFp < 0.0002 or < 0.1%. The Yp from 
this code has been shown to agree within this margin 
with somewhat independent calculations [l^. We note 
that its results are AYp ~ 0.0004 higher than those re- 
ported in Ref. [lB| for reasons that are unclear - each 
calculation should have a higher precision than this. For 
comparison, the uncertainty arising from the currently 
recommended t„ is ~ 0.0003 (but see below for its likely 
underestimation) . 

BBN codes integrate a set of coupled rate equations 
to evolve light-nuclide abundances in the expanding and 
cooling universe. The rate coefficients in those equations 
are in turn based on nuclear cross sections, mainly in- 
ferred from experiment. Several improvements in the nu- 
clear in put s have occurred in the decade since the study 
of Ref. hj, so we have updated the input database used 
there accordingly. We now summarize those changes with 
an emphasis on the D and ^He yields that are useful 
for constraining A'off. In addition to the reactions af- 
fecting Yp and D/H, the input database for the reaction 
'^IIe(a, 7)^Be has also been brought up to date by includ- 
ing recent data examined in Refs. [l^ [1^; this affects 
only the curves shown in Fig. [5] and almost none of the 
discussion. 



There are very few cross section data for neutron cap- 
ture on protons at energies relevant for BBN, though 
the situation has improved in recent years [l^] , and their 
precision is generally rather low. For this reason, BBN 
evaluations in the 1990s [13, HH] used theoretical cross 
sections from the ENDF/B-VI and earlier databases [l^ 
rather than fitting the data. In general, theoretical calcu- 
lations more precise than the data at BBN energies [26l - 
[28j were possible, but some calculations had only been 
carried out for thermal neutrons, and uncertainties were 
poorly quantified. 

The BBN calculations of Refs. [13, ^ applied to the 
ENDF cross section a conservative 5% error in overall 
normalization that had a somewhat murky origin in the 
literature. It has subsequently been shown using pionless 
effective field theory (EFT) that the low-energy cross sec- 
tion can be computed from a small number of input quan- 
tities (most at least implicitly used in the ENDF cross 
section) and shown to have an error of < 1% at BBN 
energies [1^ i30|] . We have accordingly used the EFT re- 
sults [1^ to compute the rate of p{n, j)d, assigned it an 
overall normalization error of 1%, and incorporated the 
resuhs into the BBN code in place of the ENDF/B-VI 
rate. Since the EFT rate agrees with the old ENDF rate, 
this results in a much smaller error on D/H but does not 
shift the most-likely yield. 



2. The reactions d(d, n)^He and d(d,p)^H 



The investigations in Ref. [17| resulted in a call for 
improved data on deuteron-deuteron reactions at 100 to 
a few hundred keV (center-of-mass) , where there were no 
modern data with well-quantified errors. The response 
was the work of Ref. |31j . which produced data with 
a precision of ~ 2%. We added these results to our 
rate database along with lower-energy data ^32] that were 
missed in the literature search of Ref. [l3| . The addition 
of these data reduces the error on D/H, particularly at 
low values of flbh'^ where d + d reactions formerly domi- 
nated the error budget [l3|. The new data shift D/H by 
less than 1% relative to previous results [itI. [33|. 



3. The reaction d{p, 7)^He 

It was found in Ref. [l7| and by other authors [H, 113, 
im that the only reactions important for the uncertainty 
on the D/H yield were the three reactions discussed above 
and proton capture on deuterium. Even before the above 
updates to pin, j)d and d + d rates, the d{p, 7)'^He rate 
dominated the error on D/H at the CMB-inferred ftbh^. 

The effect of using measured d{p,'-f)^Re cross sections 
in BBN calculations has not changed since our previous 
study Ref. jTJl . There has been one published low-energy 
cross section measurement since then 13611, and we have 
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updated our database to include it. As the new data 
are well below most of the BBN energy window shown in 
Fig. [H and they agree with previous data there, they do 
not affect our results discernably. This calculation will 
be referred to as "Empirical dp7" below. 

Fig. [T] shows the scarcity of laboratory data for this 
reaction at BBN energies that was cited in Ref. [17] as 
a cause for concern in both rate determination and er- 
ror estimation. Cross section fitting for charged-particle 
reactions is generally done with s-wave Coulomb barrier 
penetration divided out, in terms of the 5-factor 

S{E) = Ee^''^'^^°"'/''a{E), 

where a is the cross section, E the energy, Zi the nu- 
clear charges, a the fine-structure constant, c the speed 
of light, and v the relative velocity. Except for one rate 
evaluation using the i?- matrix formalism [s^], all BBN 
work to date has fitted S{E)eithej: as a polynomial with 
poor physical motivation flaL [igI . [25j or as an expansion 



in linearly- independent, locally-supported B-splines [17 1 
in which S{E) at a given energy is determined by nearby 
data. There is only one modern data set for d{p, 7)'^He 
at the most important energy range of 50 - 500 keV [11] , 
so only the i?-matrix approach allows independent mea- 
surements to average against it in a physically motivated 
way. 

Because ab initio nuclear theory is now quite accu- 
rate for many bound-state and continuum properties of 
three-body systems 



41[, a purely empirical approach 



to 7)^He is no longer required; the situation with 
d(p, 7)'^He is now similar to that with p{n, 7)(i. In par- 
ticular, the Pisa group has computed the cross section 
of d{p, 7)'^He using accurate two- and three-nucleon po- 
tentials, corresponding two- and three-body corrections 
to the electromagnetic currents, and accurate methods of 
solution [H, 113] . The results are in very good agreement 
with data above and below the BBN energy range (par- 
ticularly the new low-energy data of Ref. ^36] ) , not only 
in total cross section but also in angular-distribution and 
polarization observables. However, the total cross section 
is well above most of the sparse data at BBN energies. 

The energy dependence of the cross section is deter- 
mined by competition between s- and p-wave capture, 
which is probed by the polarization observables. Since 
both they and the low-energy total cross sections sup- 
port the model, we suspect that the data of Ref. [3^ 
may have a systematically low normalization by about 
20% (compared with their stated systematic error of 9%) 
and that the theoretical curve provides a better guide to 
the actual cross section. 

For these reasons, we adopt the ab initio cross section 
[39I l40j for our BBN calculations. The precision of the 
calculation is likely a few percent at worst, but it is dif- 
ficult to determine. Errors could in principle be prop- 
agated through the constraints on "model-dependent" 
two-body currents, but these are themselves small cor- 
rections to the main two-body corrections that make up 
~ 15% of the cross section. The remaining errors arise 
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FIG. 1: (Color online) Upper Panel: S'-factor data for the 
reaction d(p, 7)^He are shown, along with the best-fit curve 
and 2a intervals based purely on those data [1^ (dash-dot 
and dash-dot-dot-dot curves) and on the ab initio calcula- 
tion [39I . |40| | (solid) . All errors are shown as 2a intervals for 
comparability with Fig. 3 of Ref. [13] • Lower Panels: Sen- 
sitivity functions [l3] showing the relative importance of this 
reaction at varying energy for D/H ((72) and Li/H ((77) yields, 
evaluated at Qth'^ ~ 0.019. Convolving the S-factor error 
band with the sensitivity function gives the total uncertainty 
arising from this reaction. The strongest sensitivity to the 
dlpjjY^He rate occurs at the largest \gi\, in the 50-500 keV 
range. 



from the numerical method of computation (probably 
negligible) and the nuclear interaction model. The meth- 
ods that have been used to construct the modern nucleon- 
nucleon potentials (~ 40 parameters fitted to thousands 
of data points from many experiments) have not been 
adapted to facilitate error propagation. The modern po- 
tentials do reproduce many properties of two- and three- 
nucleon systems [2l,|4l|,|43 at the percent level, and they 
have also been applied very successfully to larger systems 
(e.g.. Hi). 

Because assigning a quantitative error to the computed 
cross section is difficult, we would like to be conserva- 
tive. The case of p{n, j)d is again analogous: through 
the 1990s, errors of 5% were assigned to an evaluation 
that later proved accurate to something more like 1% 
fvA [25I . The review of Ref. [23^ presents a polynomial 
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fit to the d{p, 7)^He S'-factor and a recommended cross 
section primarily for use well below the BBN range. This 
fit's normalization is dominated by the precise and copi- 
ous low-energy data, which cause the recommended 5'(0) 
to have an error of 7%. The ab initio cross section is in es- 
sentially perfect agreement with the low-energy data, so 
we adopt 7% as the precision to which it has been tested 
and apply this as an overall error on its normalization. 
This error is handled in our BBN Monte Carlo procedure 
just like that for the p{n, j)d cross section. Ref. [l3| in 
effect propagated a 9% error bar from the nuclear data, 
so adopting 7% does not drastically change the (domi- 
nant) contribution of d{p, 7)'^He to the yield errors. We 
use these results rather than the "Empirical dp7" calcu- 
lation described above in our likelihood analysis. 

The larger ab initio cross section for d{p, 7)^He at BBN 
energies results in more destruction of deuterium at late 
times than was found previously [13, SI] . At the WM AP 
baryon density [9j], this amounts to about a 5.6% reduc- 
tion of D/H; the D/H-inferred Ob/i^ is shifted downward 
by roughly half a standard deviation. In light of the 
severity of this shift, independent cross section measure- 
ments at 50-500 keV would be quite useful. The shift in 
the d{p, 7)'^He rate also increases the ''Li yield by about 
10% through its effect on the availability of '^He at late 
times, worsening slightly the primordial ^Li problem re- 
viewed in Ref. 1451 . 



4- The neutron Itfettme t„ 

The mean lifetime t„ of the neutron is used in the BBN 
calculation to normalize all of the weak-interaction rates 
that interconvert neutrons and protons at early times. It 
may be conveniently thought of as determining the weak 
coupling constant Gp, and its uncertainties propagate 
through the weak rates into the abundance yields. 

The work of Refs. [13, HH used a neutron mean lifetime 
of Tn = 885.4 ± 2.0 s, which was the most recent experi- 
mental result at the time. This result was subsequently 
included in the world average of 885.7 ± 0.8 s recom- 
mended for several years by the Particle Data Group [i^ . 
More recently, conflicting lifetimes of 878.5 ± 0.7 ± 0.3 s 
and 880.7 ± 1.3 ± 1.2 s have been reported respectively 
by Serebrov ct al. [13, El] and Pichlmaier et al. [4§]. The 
Particle Data Group now recommends a world average 
that includes confiicting values, 881.5 ± 1.5 s [50], with 
errors that have been inflated to reflect the discrepancy. 

We performed our calculations using the currently rec- 
ommended Tn with its recommended error, and these 
form the basis of our reported results. We also performed 
BBN calculations with the formerly recommended t„ = 
885.7±0.8 s and with t„ = 878.5±1 s from Serebrov. The 
difference between the two amounts to AYp ~ 0.0015, 
equivalent to AA'eff — 0.12. The effect on D/H is a shift 
of 0.4% at fifc/i^ = 0.0225. For both yields, the differ- 
ence is much smaller than the other uncertainties in our 
analysis. 



B. The NcB dependences of D/H and Yp 

We computed BBN yields, including their errors and 
correlation coefficients, on a grid covering flbh^ from 
0.005 to 0.1 and N^s from 0.046 to 10.046. This in- 
cludes the entire range sampled in our Markov chain 
Monte Carlo (MCMC) analysis below, allowing consid- 
erable margins at the edge of the grid. We show all BBN 
yields as functions of over this grid in the left panels 
of Fig. [2] Separate curves are shown, with widths indicat- 
ing la errors, at integer intervals in N^ff. This indicates 
the nature and strength of the sensitivity to A'cff. The 
dependence of Yp and D/H on Ncs is apparent; that of 
Yp is much larger, relative to its yield uncertainties, than 
that of D/H. At the band indicating the current WMAP 
value of rife/i^, D/H changes by about 10% per unit Ncs- 

The right side of Fig. [2]shows the dependences of BBN 
yields on Nes and on the assumed rates of d{p,j)^He 
and of the weak processes scaled by r„. Each yield was 
computed at flbh^ = 0.02260(54), inferred below from 
a cosmological analysis of CMB fluctuations marginal- 
ized over Yp and iVeS, and the error on i}},h^ was added 
in quadrature to the nuclear uncertainties. The diago- 
nal curves are thus essentially slices along the "WMAP" 
band in the left panels of Fig. [2l The effects on D/H and 
^Li/H of switching d{p,^)^tle rates are shown as black 
curves with no error band in the right side of Fig. [21 
that on D/H corresponds roughly to AA'eff = 0.42. Sim- 
ilarly, the effects of assuming either the pre- 20 11 or the 
Serebrov values of r„ are shown as black curves without 
error bands next to the Yp yield curve. To a very good 
approximation, d{p, 7)^He does not affect Yp and t„ does 
not affect D/H. 

The yields Yp and D/H measure slightly different 
things At the start of standard BBN the popula- 
tion of relativistic particles consists of photons, electrons, 
positrons, and three species of neutrinos, all at the same 
temperature. At a given temperature, the resulting ex- 
pansion rate of the universe is 2.3 times that from pho- 
tons alone. Weak freezeout begins at this time, setting 
the neutron/proton ratio that will eventually determine 
Yp. An additional neutrino species speeds the expansion 
rate by ~ 40.3%, so that weak freezeout occurs at higher 
temperature. Higher temperature at equilibrium implies 
more neutrons, so faster expansion yields more ^He. 

In contrast, the abundance of deuterium is determined 
by c? -I- d and d + p reactions toward the end of BBN, 
when the photon temperature is well below the electron 
rest mass. At these times, there are essentially no elec- 
trons or positrons, and the energy they previously carried 
has been transferred to the photons, which are conse- 
quently at higher temperature than the neutrinos. The 
expansion rate at this time in standard BBN is ~ 1.7 
times that of photons alone, so that adding a neutrino 
species at the same temperature as the others causes a 
~ 36.7% speed-up. Faster expansion implies less time 
available for deuterium burning and thus higher D/H. In 
the standard model, deuterium burning and the recom- 



5 




FIG. 2: (Color online) Left panels: Yields of light nuclides as functions of the baryon density, Q,bh^ ■ Each band indicates yields 
for a single value (integer plus 0.046) of A^eff and its thickness the Icr nuclear uncertainty on those yields. Bands are shaded 
lighter for smaller values of NeS and darker for larger A'eff . Also shown are horizontal bands indicating observational constraints 
on abundances [5ll - [55| and a vertical band indicating the WMAP7 la interval for Q,hh? ■ Right panels: The same yields shown 
as functions of A'eff at the CMB-inferred value of Qbh^, essentially slices along the band labeled "WMAP" in the left panels. 
Salmon-colored bands indicate yields with our adopted rates, and their widths indicate a quadrature sum of nuclear errors and 
the error on Q,i,h^ . In the Yp graph, blue la and 2a contours indicate constraints from WMAP7-I-SPT alone, while black curves 
paralleling the adopted yields indicate the effects of using the old Particle Data Group [i^ (upper) and Serebrov [i^ (lower) 
values of r„ (no errors shown). In the D/H and ^Li/H graphs, the thin solid curves indicate the result of using the empirical 
rate for d(p, 7)^He, again without showing the errors. Horizontal bands indicate observational constraints as discussed in the 
text. 



bination era share the same populations of relativistic 
species: photons plus three neutrino species at (4/11)^/'^ 
of the photon temperature. 



C. Observed primordial abundances 



Constraints based on BBN require observed abun- 
dances. The inference of each primordial abundance is 
a highly technical subject unto itself with its own exten- 
sive literature. Here we provide brief descriptions of the 
experimental constraints and explain the specific values 
that we choose. The interested reader is referred to the 
cited literature for more information, and to the reviews 
inRefs. i, IH, [HO, IHB-IEi . 



'He 



The primordial ^He mass fraction Yp is inferred from 
emission-line spectra of Hii regions in metal-poor blue 
compact galaxies. Lines of H and He are observed and a 
He/H ratio inferred using a radiative transfer model; the 
result is customarily expressed in terms of Y-p. Because of 
the weak dependence of Yp on VL^h? evident in Fig. [21 the 
desired precision of measurement is in the percent range. 
Large data sets have been assembled, and the errors on 
Y-p are dominated by systematic effects in the analysis, 
not statistics. 

A recent revision of the atomic data used to infer Y-p 
[59| resulted in a systematic upward shift of all data sets 
(but not all individual data). Three groups have pub- 
lished Yp values based on the new emissivities, and we 
mention only the most recent publication from each. Pe- 
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imbert et al. [511 find 0.2477 ± 0.0029 from 5 Hii regions 
using radiative-transfer models tailored to each object 
and taking into account temperature fluctuations within 
each Hii region. Applying approaches in which physical 
parameters for each Hii region are determined in a si- 
multaneous fit and errors are computed by Monte Carlo, 
Izotov & Thuan find 0.2565±0.0010 (stat)±0.0050 (syst) 
from 93 Hii regions j54i] while Aver et al. find 0.2573 ± 
0.0033 [ill for nine selected low-metallicity objects from 
the Izotov & Thuan data set. 

The quoted Yp is generally an extrapolation to zero 
metallicity Z of the helium mass fractions Y observed in 
objects of (small) nonzero Z. Peimbert et al. applied a 
dY/dZ slope calibrated elsewhere; Izotov & Thuan deter- 
mined dY/dZ internally to their data by linear regression; 
the Aver et al. result quoted above is the mean of their 
nine selected objects, with zero slope assumed. They also 
derived from linear regression with no prior on dY/dZ a 
result of 0.2609±0.0117, but it appears that their narrow 
metallicity range does not allow a slope to be determined. 
The upper end of this error interval is close to the ini- 
tial Pop I composition of the Sun, Fproto = 0.278 ± 0.008 
[60j . which must be greater than Y-p since stars are net 
creators of ''He. We repeat our analysis for each of these 
three input Yp constraints. 



2. 



'-H 



Primordial deuterium is observed in gas of very low 
metallicity along our lines of sight to distant quasars. 
Quasar light is absorbed by atomic transitions of neutral 
hydrogen, and deuterium is observed as additional ab- 
sorption with the appropriate isotope shift and doppler 
width to be deuterium in the same physical location. 
Quasar spectra that are suitable for inference of D/H 
are rare, because the quasar continuum emission must 
be well-characterized, there must be no other absorption 
nearby in redshift/velocity space, and the absorber must 
have sufficient hydrogen column density to produce mea- 
surable absorption from its small amount of deuterium. 
To date there are eight extragalactic systems in which 
D/H is generally agreed to have been measured. These 
are described in Refs. [1^ |6l| (and references therein), 
where the average log(D/H) = —4.556 ± 0.034 is found. 
As discussed in Refs. [53, H^, 113, it is not clear that all 
eight systems have consistent inferred D/H. We adopt 
this average as the observed primordial D/H. 



3. ^He and Li 

The other two nuclides predicted in measurable quan- 
tities by BBN are ^He and ^Li. ^He has only been ob- 
served in our galaxy (in Hii regions) and only at rel- 
atively high metallicity 6J|. Moreover, its post-BBN 
history is complicated, as it seems to be produced co- 
piously in some stars but mostly burned up in oth- 



ers [63. Although its observed abundance is generally 
consistent with BBN, further inferences from ^He seem 
unwarranted. The primordial ^Li/H has been identi- 
fied with the nearly-uniform abundance found in atmo- 
spheres of low-metallicity main-sequence turnoff stars of 
the Galactic halo [1^, H^, IH, [13 . As shown in the left 
side of Fig. [51 this abundance is near the minimum pro- 
ducible by standard BBN at any ilj, and a factor of 3.5 
below that predicted by the WMAP fif,. Whether this 
is a consequence of post-BBN processing (most likely in 
the observed stars themselves) or of something missing 
in standard BBN is unclear Since we are consid- 

ering only standard BBN apart from variation of Ncs, 
the important fact here is that even at N^s = 10, BBN 
overproduces ''Li by a factor of 2 relative to halo stars. 
Thus, it is an implicit assumption of our model that the 
observed Li/H is not primordial, and we do not include 
it in our constraints. 



III. BBN AND CMB LIKELIHOOD FUNCTIONS 

The cosmic microwave background is affected by extra 
radiation in two ways On large scales the CMB 

is sensitive to the gravitational effect of the perturba- 
tions in the extra radiation, while on smaller scales the 
increased expansion rate leads to an increased amount of 
diffusion damping. However, changing the helium abun- 
dance also affects diffusion damping, as Yp affects the 
number of electrons per unit baryon mass. The recent 
results of Keisler et al. [ll| used CMB power spectra 
from the WMAP seven-year data release (WMAP7) and 
the South Pole Telescope (SPT) to find A^eff = 3.85±0.62 
when the relation between Yp and A^eff from BBN theory 
was enforced, ignoring uncertainties in the BBN rates. 
Without an imposed BBN prior on Yp the precision of 
the A'ofr measurement is degraded to A'cfr = 3.4 ± 1.0. 
Following Keisler et al. [il[, we use bandpowers from 
WMAP7 and SPT, which together cover the multipole 
range £ = 2 - 3000. 

To combine BBN constraints with CMB measure- 
ments, we adopt a Markov Chain Monte Carlo (MCMC) 
approach, using a modified version of the publicly avail- 
able package CosmoMC [3 • MCMC methods work by 
sampling the parameter space at random positions and 
either accepting or rejecting each sample based on the 
likelihood of that particular point in the parameter space. 
The most-likely values and confidence intervals for the 
parameters are identified from the distribution of sam- 
ples. 

For CMB measurements, we assume an 8-parameter 
spatially flat ACDM cosmology, allowing as parameters 
the cold dark matter density {flch'^), the baryon den- 
sity (ilfcft.^), the angular size of the sound horizon {0s), 
the optical depth to Thomson scattering in the reionized 
universe (r), the amplitude (Ag, using as a parameter 
In As) and spectral index of the potential fluctuations 
(jig), the helium abundance {Yp) and the effective num- 
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ber of neutrino species (A^cff)- As additional parameters, 
we allowed for a contribution to the SPT power from 
the Sunyaev-Zeldovich effect, clustered galaxies, and shot 
noise from galaxies, using priors from ShirokofF et al. 69], 
as described in Keisler et al. We assume broad 

uniform priors that extend well beyond the WMAP7- 
allowed range on all cosmological parameters and explore 
< < 1 and 1.046 < iVcff < 8.046. 

The BBN likelihood only depends on a small subset 
of these parameters: flbh"^, N^s, and Yp. For a given 
assumed pair of fl^h^ and N^s, a probability distribution 
based on nuclear uncertainties can be calculated for Yp 
and D/H from BBN theory. This can be multiplied by 
the observed joint probability distribution describing Yp 
and D/H observations, and the resulting distribution can 
then be integrated over D/H to yield a final P{Yp): 

P{Yp, D/H, N,s, nt,h^\Yp^^^,D/Hobs) cc 
PiYp^bs^D/Hobs\Yp,D/H) 

xP{Yp,D/H\N,s,^bh^) (1) 

As written, the two factors on the right side of the equa- 
tion are the observed abundances and the theory predic- 
tion for the abundances, each of which is a 2-D probabil- 
ity distribution. 

The D/H ratio does not affect CMB observables, so we 
marginalize over it to obtain the final likelihood: 

P{Yp,N,s,nbh^)oc J PiYp,D/H,N,if,nbh^)d{D/H) 

(2) 

In general, this procedure could be time-consuming in a 
MCMC, where likelihood evaluations may be done hun- 
dreds of thousands of times. We approximate the ob- 
served and theory abundances as 2D Gaussians, with the 
observed abundances uncorrelated and the (small) the- 
ory covariances extracted from the Monte Carlo BBN 
calculations described above. In this case, Eq. ([T]) is 
simply the product of two Gaussians, the result of which 
is another Gaussian distribution, and the marginaliza- 
tion is trivial. Therefore, the likelihood evaluation and 
marginalization over the true D/H are reduced to an al- 
gebraic equation involving the means and covariances of 
the observed and theory abundances. 

In more detail, following the notation of Ahrendt 70], 
if the probability distributions of the theoretical and ob- 
servational abundances are given by multivariate normal- 
ized normal distributions N{pLth, "^th) and N{fIobs, '^obs), 
with means denoted by fl and covariance matrices E, the 
product of these two distributions is 

N{flth, ^th)N{flobs, ^obs) = {^l^corabi ^comb ) (3) 

where the combined covariance matrix is 



'^comb — 



-1 
th 



and the combined mean is 

f^comb 



^comby^th 

fith + ^obsi^obs 



(4) 



(5) 



The prefactor is given by 

_ 1 
" |27r(Ei„ + I]ob,)|-i/2 
1 



(6) 



exp 



{p-th — fiobs)^{^th + 5]o6s) ^{fith — j^obs) 



Marginalizing over D/H is then done by an analytic inte- 
gral which leads to the 2-D normal distribution in Equa- 
tion (3) being simply replaced by a 1-D normal distribu- 
tion involving only the mean and covariance of Yp . 

In summary, to include full BBN information within 
an MCMC, one must simply supply observed He and D 
abundances and a covariance matrix for the observations 
(presumably always diagonal) , and a grid of BBN predic- 
tions as functions of ftbh^ and A^cff for He and D abun- 
dances with the theoretical uncertainty expressed as a 
covariance matrix. While not pursued here, the general- 
ization to include other abundances (such as "^He or ^Li) 
is straightforward in this framework. 



IV. BBN AND CMB CONSTRAINTS ON N^s 

Before combining BBN and CMB constraints, it is im- 
portant to first check for consistency. CMB-only con- 
straints on Nos marginalized over cosmological parame- 
ters and Yp are weak but still useful, as can be seen in 
Fig. m and Table U (and more fully in Fig. 12 of Ref. 
[ll|). Furthermore, we can use measurements of helium 
abundances to improve the precision of the CMB mea- 
surement. This is true even with weak constraints on Yp, 
which can be provided by the protosolar abundance as 
an upper limit. 

A lower limit on Yp has less impact but also helps. 
It has long been recognized that observed helium abun- 
dances everywhere are too large to be attributed to stars 
[tH . Both observations and models of chemical evolution 
find AY/AZ in the range 1 to 5 {72,] at low metallicities 
and smaller slopes at higher metallicities where type la 
supernovae contribute to Z. A solar metallicity of 0.014 
[73j . would then imply that stars have contributed at 
most Ay ~ 0.07. Since Yp^^to = 0.278 ± 0.008 is the 
protosolar abundance, this implies Yp > 0.22, consistent 
with all observations over a long period of time (cf. Table 
1 of Ref. [zll and the accompanying discussion). 

Applying the weak prior that 0.22 < Yp < l^proto, the 
precision of the CMB measurement of TVcff is improved to 
3.87 ± 0.81. This is only slightly weaker than constraints 
on Ncs using CMB data combined with the substantially 
more precise measurements of Yp from nearby galaxies, 
which are shown as the first column in Tabic HI 

Separately, the more precise determinations of Yp 
can be combined with D/H measurements and BBN 
theory to get a BBN-only measurement of A^off (and 
a simultaneously-fit flbh'^ that agrees with the CMB), 
shown in the last column in Table [H The constraints are 
formally quite strong, at AA'cft ^ 0.25, but results using 
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different Yp determinations inherit the sometimes > 2a 
disagreement of the underlying Yp . 

A direct comparison of the first and last data columns 
provides a non-trivial test of cosmology. For the two 
higher Y-p determinations, the separate iVcff values from 
BBN and the CMB are in excellent agreement, suggest- 
ing that the effective number of neutrino species was the 
same in the first few minutes of the early universe as it 
was at ~ 300 000 years. The lower Y-p measurement of 
Peimbert is in mild tension with this interpretation, with 
only a very slight preference for an extra neutrino species 
at late times but better agreement with the expectation 
that iVcff = 3.046. 

Given the apparent consistency, it is reasonable to 
combine CMB and BBN constraints. Simply requir- 
ing that the relation between Yp and NcS be consistent 
with standard BBN theory leads to an improvement in 
the CMB measurement, and further combining this con- 
straint with observed light element abundances leads to 
the very tight constraints shown in the second and third 
data columns. Again, the systematically different Yp 
determinations lead to challenges in interpretation: the 
combined analysis yields either a clear detection of an 
excess of effective neutrinos over the standard value of 
3.046 (at 2.2 — 3.4(t) or excellent agreement (within lu). 

Given the disagreement between 1^ determinations, it 
is interesting that a constraint on A'eff can be obtained 
using CMB measurements, BBN theory, and only the 
D/H abundance, as can be seen in Figure[3] Using no Yp 
measurements at all, there is a hint of excess radiation (at 
just below 2<t). This is interesting, as the current mea- 
surement of D/H is at the 10% level. It is challenging to 
find systems that are well-suited for this measurement, 
but it is not inconceivable that the D/H abundance mea- 
surement could be sustantially improved. D/H has long 
been viewed mainly as a probe of fi&/i^, but it may have 
a future as a probe of the expansion rate encoded by A'eff 
with independent systematics, separate from Yp. 



V. CONCLUSION 

In this paper, we have presented a comparison of de- 
terminations of the effective number of neutrino species, 
TVeff, from recent CMB measurements (WMAP7 and 
SPT) and BBN calculations, applying observed abun- 
dances of helium and deuterium in several ways. Our 
results are similar to those of other analyses in the recent 
literature ItsI-Itt'I , though comparison is difficult be- 
cause of the use of different cosmological data and model 
spaces and the common practice of setting Yp — 0.24, 
which seems on the basis of Fig. [2] to bias N^s a bit 
high. Like these and previous analyses, ours generally 
favor iVefi > 3. The most directly comparable result in 
the literature combines BBN and CMB data with con- 
straints on large-scale structure and the Hubble constant 
to find 95% confidence limits of iVcff = 3.90l:|^j^ 
This is similar to our results. 




FIG. 3: (Color online) Constraints on A'eff from just CMB 
and D/H measurements, shown as marginalized likelihoods. 
The black solid curve is CMB alone, the black dotted curve 
shows CMB with BBN constraints on the relation between 
Vp and A'eff , and the red dashed curve also adds the observed 
abundance of D/ H . 



The BBN calculation used here includes a new estimate 
of a key cross section, (i(p, 7)'^He, based on detailed new 
theoretical calculations, which shifts deuterium yields by 
~ 6% and lithium yields by ~ 10%. This new cross 
section shifts Vlbh^ or Yp inferred from D/H and slightly 
aggravates the large discrepancy between standard BBN 
and observed ''Li abundances. 

The BBN calculations and observed abundances have 
been embedded in a module that can be used in a widely 
used MCMC code (CosmoMC) and combined with CMB 
measurements or other cosmological data. The BBN like- 
lihood calculation includes both theoretical and observa- 
tional uncertainties, including the weak correlation be- 
tween computed D/H and Yp. 

We found that there is good agreement between BBN- 
based and CMB-based determinations of A'eff. These 
tests probe two very different epochs of the universe (a 
few minutes versus a few hundred thousand years), pro- 
viding evidence that the same cosmological model de- 
scribes both epochs. However, both probes show hints 
of new physics in that the preferred number of effective 
neutrino species is above the standard model expectation. 
This is especially true in the case of the two higher Yp 
determinations used (based on overlapping data), with 
as much as a 3.4(t excess. 

Setting aside the helium constraints and combining the 
CMB with BBN theory and the observed D/H, the ex- 
cess is 2(7, at A'eff = 3.9 ± 0.44. Improvements in the 
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measurement of the D /H abundance would be extremely 
valuable as both a cosmological probe and an external 
check on Yp determinations. 

The results here underscore the complementarity of 
CMB measurements, BBN theory, and light element 
abundance measurements. CMB measurements allow 
both an increased precision in BBN predictions (through 
a robust measurement of the baryon density) and a strong 
test for BBN (through comparisons with A'eff and Yp de- 
terminations). This is a consequence of the strong link- 
ages in the standard cosmology between light element 
yields in the early universe and observable signatures at 
recombination, which connect the underlying cosmolog- 
ical model, the acoustic oscillations and Silk damping 
observed in the CMB, and the abundances observed in 
low-metallicity systems. Precise measurements of all of 
these interlocking observables provide strong tests of the 
model and its parameters. 
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3.87 ±0.24 
3.83 ± 0.35 
3.22 ±0.19 


3.85 ± 0.26 

3.82 ±0.45 
3.13 ±0.21 


CMB only 
CMB+0.22 < Yp < Yproto 
CMB+BBN 
CMB+BBN+D/H 


3.4 ±1.0 
3.87 ± 0.81 

3.89 ± 0.60 

3.90 ± 0.44 



TABLE I: Constraints on iVefi from various combinations of data sets. The first three rows of results use three different precise 

determinations of Yp, while the corresponding columns show different combinations of these Yp measurements with CMB data, 
BEN theory (including errors on reaction rates), and D/H abundance measurements. The lower set of results do not use these 
precise Yp measurements. 



